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ABSTRACT 

A photometric sample of ~7100 V < 25.3 Lyman break galaxies (LBGs) has been selected by 
combining Subaru/Suprime-Cam BVRci'z' optical data with deep G^LiJX/ A^f/l^ imaging of the 
Subaru Deep Field. Follow-up spectroscopy confirmed 24 LBGs at 1.5 z < 2.7. Among the optical 
spectra, 12 have Lya emission with rest-frame equivalent widths of « 5 — 60A. The success rate for 
identifying LBGs as iVf/F-dropouts at 1.5 < z < 2.7 is 86%. The rest-frame UV (1700A) luminosity 
function (LF) is constructed from the photometric sample with corrections for stellar contamination 
and z < 1.5 interlopers (lower limits). The LF is 1.7 ± 0.1 (1.4 ± 0.1 with a hard upper limit on 
stellar contamination) times higher than those of z '--^ 2 BXs and z ~ 3 LBGs. Three explanations 
were considered, and it is argued that significantly underestimating low-z contamination or effective 
comoving volume is unlikely: the former would be inconsistent with the spectroscopic sample at 93% 
confidence, and the second explanation would not resolve the discrepancy. The third scenario is that 
different photometric selection of the samples yields non-identical galaxy populations, such that some 
BX galaxies are LBGs and vice versa. This argument is supported by a higher surface density of LBGs 
at all magnitudes while the redshift distribution of the two populations is nearly identical. This study, 
when combined with other star-formation rate (SFR) density UV measurements from LBG surveys, 
indicates that there is a rise in the SFR density: a factor of 3 — 6 (3 — 10) increase from z ~ 5 (z ~ 6) 
to z ^ 2, followed by a decrease to z ^ 0. This result, along with past sub-mm studies that find a 
peak at z ~ 2 in their redshift distribution, suggest that z 2 is the epoch of peak star-formation. 
Additional spectroscopy is required to characterize the complete shape of the z ~ 2 LBG UV LF via 
measurements of AGN, stellar, and low-z contamination and accurate distances. 
Subject headings: galaxies: photometry — galaxies: high redshift — galaxies: luminosity function — 
galaxies: evolution 



1. INTRODUCTION 

Over the past decade, the n umber of Lyman break 
galaxies (LBGs; for a review, see lGiavaliscci|[200^ identi- 
fied at z ~ 3 — 6 has grown rapidly fr om deep, wide- 
field optical imaging surveys (e.g., Steidel et al.l 119991 : 
iBouwens et"all l2006t lYoshida et all l2006f l. Follow-up 
spectroscopy on large telescopes has shown that this 
method (called the Lyman break technique or the "drop- 
out" method) is efficient at identifying high-z star- 
forming galaxies. Furthermore, these studies have mea- 
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sured the cosmic star-formation history (SFH) at z > 3, 
which is key for understanding galaxy evolution. It indi- 
cates that the star-formation rate (SFR) density is 10 or 
more times higher in the past than at z ~ 0. 

Extending the Lyman break technique to z < 3 re- 
quires deep, wide-field UV imaging from space, which is 
difficult. In addition, [O II] (the bluest optical nebular 
emission line) is redshifted into the near-infrared (NIR) 
for z > 1.5 where high background and lower sensitiv- 
ity limit surveys to small samples (e.g., Malkan et al] 
1996; 'Moorwo od et all 120001 : Ivan der Werf et all l20dol : 
Erb et al. 2003). The combination of these observational 
limitations has made it difficult to probe z « 1.5 — 2.5. 

One s olution to the prob l em is the 'BX' method devel- 
oped bv lAdelberger et all (|2004D . This technique iden- 
tifies blue galaxies that are detected in U, but show a 
moderately red U — G color when the Lyman continuum 
break begins to enter into the ?7-band at z '-^ 2. 

Other methods have used NIR imaging to identify 
galaxies at z — 1 — 3 via the Balmer/4000A break. 
For example, selection of objects with J — K > 2.3 
Vega) has yielded "dista nt red galaxies" at z ~ 2 — 3 
van Dokkum et al] |2004| ) . and the ^BzK method has 
found passive and star -forming (dusty and less dusty) 
galax ies at z « 1.5—2.5 (|Daddi et al.ll2004tlHavashi et al.l 
|2007( ). The completeness of these methods is not as 
well understood as UV-selected techniques, since limited 
spectra have been obtained. 

In this paper, the Lyman break technique is extended 
down to z '--^ 1.8 with wide-field, deep A^[/V^ imaging of 
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the Subaru Deep Fi eld (SDF) with the Galaxy Evolution 
Explorer (G^L£^X: [Martin et al.lf2005l ). This survey has 
the advantage of sampling a large contiguous area, which 
allows for large scale structure studies (to be discussed 
in future work), an accurate measurement of a large por- 
tion of the luminosity function, and determining if the 
SFH peaks at z ~ 2. 

In § [21 the photometric and spectroscopic observa- 
tions are described. Section [3| presents the color selec- 
tion criteria to produce a photometric sample of NUV- 
dropouts, which are objects undetected or very faint in 
the NUV, but present in the optical. The removal of 
foreground stars and \ow-z galaxy contaminants, and 
the sample completeness are discussed in § [H In § [SJ 
the observed UV luminosity function (LF) is constructed 
from ~7100 A^f/V^-dropouts in the SDF, and the comov- 
ing star-formation rate (SFR) density aX z — 1.8 — 2.8 is 
determined. Comparisons of these results with previous 
surveys are described in § [Bl and a discussion is provided 
in § [71 The appendix includes a description of objects 
with unusual spectral properties. A flat cosmology with 
[r^A, r^M, ^7o] = [0.7, 0.3, 1.0] is adopted for consistency 
with recent LBG s tudies. Al l magnitudes are reported 
on the AB system (|Okelll974D . 

2. OBSERVATIONS 

This section describes the deep NUV data obtained 
(ij 12. ip . followed by the spectroscopic observations 12.21 
and [13) from Keck, Subaru, and MMT (Multiple Mir- 
ror Telescope). An objective method for obtaining red- 
shifts, cross-correlating spectra with templates, is pre- 
sented (i) l2.3p and confirms that most A^[/V-dropouts are 
at ~ 2. These spectra are later used in ^ 13.21 to define 
the final empirical selection criteria for z ~ 2 LBGs. A 
summary of the success rate for finding z ~ 2 galaxies as 
A^t/F-dropouts is included. 

2.1. GALEX/NUV Imaging of the SDF 

The SDF (jKashikawa et all I2004D . centered at 
a(J2000) = 13'^24'^38!9, (5(J2000) = -h27°29'25'.'9, is 
a deep wide-field (857.5 arcmin ) extragalactic sur- 
vey with optical da ta obtained from Suprime-Cam 
([Mivazaki et al.|[2002D . the prime-focus cam era mounted 
on the Subaru Telescope (jive et al.l I2004D . It was im- 
aged with GALEX in the NUV (1750 - 2750A) between 
2005 March 10 and 2007 May 29 (GIl-065) with a total 
integration time of 138176 seconds. A total of 37802 
objects are detected in the full NUV image down to 
a depth of «27.0 mag (3cr, 7.5" diameter aperture). 
The GALEX-SDF photometric catalog will be presented 
in future work. For now, objects undetected or faint 
{NUV > 25.5) in the NUV ure discussed. 

The NUV image did not require mosaicking to cover 
the SDF, since the GALEX field-of-view (FOV) is larger 
and the center of the SDF is located at (-H3.87', +3.72") 
from the center of the NUV image. The NUV spa- 
tial resolution (FWHM) is 5.25", and was found to 
vary by no more than 6% across the region of interest 
(jMorrissev et al.ll2007f) . 

2.2. Follow-up Spectroscopy 
2.2.1. Keck/LRIS 
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Fig. 1. — Postage stamps for some AfC/V-dropouts targeted with 
LRIS. From left to right is NUV, B, and V. Each image is 24" on 
a side and reveals that optical sources do not have a NUV coun- 
terpart. Photometric and spectroscopic information are provided 
in Table [1] 

When objects for Keck spectroscopy were selected, the 
NUV observations had accumulated 79598 seconds. Al- 
though the selection criteria and photometric catalog are 
revised later in this paper, a brief description of the origi- 
nal selection is provided, since it is the basis for the Keck 
sample. An initial A^f/F-dropout catalog (hereafter ver. 
1) of sources with NUV - B > 1.5 and B - ^ < 0.5 
was obtained. No aperture correction was applied to 
the 7.5" aperture NUV flux and the 2" aperture was 
used for optical photometry. These differ from the fi- 
nal selection discussed in § 13.21 The NUV 3a limit- 
ing magnitude for the ver. 1 catalog is 27.0 within a 
3.39" radius aperture. Postage stamps (see Figure [T]) 
were examined for follow-up targets to ensure that they 
are indeed A^[/F-dropouts. The Keck Low Resol ution 
Imaging and Spectrograph (LRIS: [Oke et al.lll995l ) was 
used to target candidate LBGs in multi-slit mode on 
2007 January 23—25. The total integration times were 
either 3400, 3600, or 4833 seconds, and 36 A^t/V^-dropouts 
were targeted within 3 slitmasks. A dichroic beam split- 
ter was used with the 600 lines mm~^ grism blazed at 
4000A and the 400 hues mm^^ grating blazed at 8500A, 
yielding blue (red) spectral coverage of 3500 - 5300A 
(6600 — 9000A) , although the coverages varied with loca- 
tion along the dispersion axis. The slits were 4" to 8" in 
length and 1" in width, yielding spectral resolution of 
«0.9A at 4300A and «1.2A at 8OOOA. 

Standard methods for reducing optical spectra were 
followed in PyRAF where an IRAF script, developed by 
K. Adelberger to reduce LRIS data, was used. When 
reducing the blue spectra, dome flat-fields were not used 
due to the known LRIS ghosting problem. Other LRIS 
users have avoided flat-fielding their blue spectra, since 
the CCD response is mostly flat (D. Stern, priv. comm). 

HgNe arc-lamps were used for wavelength calibration 
of the blue side while OH sky-lines were used for the 
red side. Typical wavelength RMS was less than 0.1 A. 
For flux calibration , long-slit spectra of BD-l-26 2606 
l|Oke fc Gunnl[T98l were obtained following the last ob- 
servation for each night. 

In the first mask, three of five alignment stars had coor- 
dinates that were randomly off by as much as 1" from the 
true coordinates. These stars were taken from the USNO 
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catalog, where as the better alignment stars were from 
the 2MASS catalog with a few tenths of an arcsecond 
offsets. This hindered accurate alignment, and resulted 
in a lower success rate of detection: the first mask had 7 
of 12 A^[/T^-dropouts that were not identified, while the 
other two masks had 2/10 and 3/14. 

2.2.2. MMT/Hectospec 

Spectra of A^f/l^-dropouts from the final photomet- 
ric catalog were obtaine d with the multifiber optical 
spectrograph Hectospec (jFabricant et all l2005f) on the 
6.5m MMT on 2008 March 13 and April 10, 11, and 
14. Compared to Keck/LRIS, MMT/Hectospec has a 
smaller collecting area and lower throughput in the blue, 
so fewer detections were anticipated. Therefore, obser- 
vations were restricted to bright (Vauto = 22.0 — 23.0) 
sources, which used 21 of 943 fibers from four configu- 
rations. Each source was observed in four, six, or seven 
20-minute exposures using the 270 mm~^ grating. This 
yielded a spectral coverage of 4000 - 9000A with 6A res- 
olution. The spectra were wavelength calibrated, and 
sky-subt racted using the stand ard Hectospec reduction 
pipeline (Fabrican t et al.l[2005| ). A more detailed discus- 
sion of these observations is deferred to a forthcoming 
paper (Ly et al. 2008, in prep.). 

2.3. Spectroscopic Identification of Sources 

The IRAF task, xcsao from the RVSAO package 
(jKurtz fc MinkI Il998l ver. 2.5.0), was used to cross- 
correlate with six UV spectral templates of LBGs. For 
cases with Lyct in emission , the composite of 811 LBGs 
from IShaplev et al.l (|2003[ ) and t he two top quartile 
bins (in Lya equivalent width) of iSteidel et al.l ()2003f ) 
were used. For sources lacking Lya emission (i.e., pure 
absorption-line system s), the spectra of M S 1512-cB58 
(hereafter 'cB58') from Pett ini et al.l (|2000| ), and the two 
lowest quartile bins of .Steidel et all (llool were used. 

When no blue features were present, the red end of 
the spectrum was examined. An object could still be 
at z > 1.5, but at a low enough redshift for Lya to 
be shortward of the spectral window. In this case, 
rest-frame NUV features, such as Fe II and Mg II, are 
available. iSavaglio et al.l (|2004D provided a composite 
rest-frame NUV spectrum of 13 star-forming galaxies at 
1.3 < 2 < 2. For objects below z « 1.5, optical features 
are available to determine redshift. Th e composite SPSS 
spectra (3500 - 7000A coverage) from lYip et all (|2004f ) 
and those provided with RVSAO (3000 — 7000A) are used 
for low-z cases. Note that in computing redshifts, sev- 
eral different initial guesses were made to determine the 
global peak of the cross-correlation. In most cases, the 
solutions converged to the same redshift when the initial 
guesses are very different. The exceptions are classified 
as 'ambiguous'. 

Where spectra had poor S/N, although a redshift was 
obtained for the source, the reliability of identification 
(as given by xcsao's i?- value) was low {R = 2 — 3). 
An objective test, which was performed to determine 
what _R- values are reliable, was to remove the Lya emis- 
sion from those spectra and templates, and then re-run 
xcsao to see what i?-values are obtained based on ab- 
sorption line features in the spectra. Among 10 cases 
(from LRIS spectroscopy), 6 were reconfirmed at a sim- 
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Fig. 2. — Cross-correlation spectra for targets that yielded R = 
2.5 — 3.2 without an emission line. From top to bottom shows 
Af;71/-dropout ID 182284, 186254, 96927, and 92942. The top two 
have ij- values of ^3.1 and are identified as LBGs while the latter 
two have R = 2.6 — 3.0 and are classified as ambiguous. The peak 
near the center of the plots represents the strongest peak in the 
cross-correlation. 

ilar redshift (Az = 2.4 x 10""* - 1.5 x lO'^)" with R- 
values of 2.30—7.07. This test indicates that a thresh- 
old of R = 2.5 is reasonable for defining whether the 
redshift of a source (lacking emission li nes) was deter- 
mined . This cut is further supported bv lKurtz fc MinkI 
(|1998D . who found that the success of determining red- 
shifts at i? = 2.5 — 3 is ~90%. However, to obtain more 
reliable redshifts, a more stricter R = 3.0 threshold is 
adopted. If a i? = 2.5 threshold is adopted, then seven 
sources with R = 2.5-3 (ID 86765, 92942, 96927, 153628, 
169090, 190498, and 190947) are re-classified as 'identi- 
fied'. These redshifts are marginally significant: a few 
to several absorption features coincide with the expected 
UV lines for the best-fit redshifts of ^2, but a few ad- 
ditional absorption lines are not evident in the low S/N 
data. Statistics presented below are provided for both 
adopted i?-value cuts. 

While some sources are classified ambiguous, it is likely 
that they could be at high-z. For example, 185177 (clas- 
sified as ambiguous) could be a LBG, since it shows a 
weak emission line at '^4500A [z ~ 2.7 if Lya) and a few 
absorption lines. This source, statistics ( ^50% success- 
ful identification for R = 2.0 — 2.5) from lKurtz fc MinkI 
(|T998), and NUV-78625 (with R = 2.3 but identified 
'by eye' to be a 2; ~ 1.6 AGN) suggest that while a 
cut is placed at R — 2.5 or R — 3.0, it could be that 
some solutions with R ~ 2.0 — 3.0 are correctly iden- 
tified. An R = 3.0 (i? = 2.5) typically corresponds to 
a peak of ~0.25 ('^0.2) in the cross correlation spectra, 
which is typically 3tT (2 — 3cr) above the RMS in the 
cross-correlation (see Figure [2]). 

2.3.1. LRIS Results 

12 (14 with R > 2.5) LBGs are found at 1.7 < z < 2.7 
out of 36 attempts. Among those, 10 show Lya in emis- 

This is lower, but still consistent with differences between 
emission and absor ption redshifts of 650 km s~^ for LBGs 
llShaplev et al.|[2003) . 
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sion, while 2 (4 with R > 2.5) are identified purely by 
UV absorption lines. Their spectra are shown in Fig- 
ures [3] and m and Table [1] summarizes their photometric 
and spectroscopic properties. Contamination was found 
from 3 stars and 5 (7 with R > 2.5) low-z galaxies (shown 
in Figure[5|), corresponding to a 60% success rate (58% if 
i? > 2.5 is adopted). Four sources showed a single emis- 
sion line, which is believed to be [O II] at z ~ 1 — 1.5, one 
source showed [O II], H/3, and [O III] at z 0.7, and two 
sources with absorption lines have i? ^ 2.5 results with 
z ~ 0.1 and ^ 0.5 (these would be "ambiguous" with 
the R > 3.0 criterion). The success of identifying z '--^ 2 
LBGs improves with different color selection criteria that 
remove most interlopers (see 13. 2p . 

Of the remaining 16 spectra (12 with R > 2.5 cut), 8 
(4 with R > 2.5 cut) were detected, but the S/N of the 
spectra was too low, and the other 8 were undetected. 
These objects were unsuccessful due to the short integra- 
tion time of about one hour and their faintness (average 
V magnitude of 24.2). 

It is worthwhile to indicate that the fraction of LRIS 
spectra with Lyg eni ission is high (83%). In comparison, 
IShaplev et all (|2003f ) reported that 68% of their z - 3 
spectroscopic sample contained Lya in emission. If the 
fraction of LBGs with Lya emission does not increase 
from z ~ 3 to z ^ 2, it would imply that 5 z ~ 2 galaxies 
would not show Lya in emission. Considering the diffi- 
culties with detecting Lya in absorption with relatively 
short integration times, the above 83% is not surprising, 
and suggests that most of the z > 1.5 ambiguous LRIS 
redshifts listed in Tabled] are correct. 

2.3.2. Hectospec Results 

Among 21 spectra, 7 objects (2 are AGNs) are iden- 
tified {R > 3.0; 9 if i? > 2.5) at z > 1.5, 2 objects are 
stars, 1 (2 with R > 2.5) is a z < 1.5 interloper, and 11 
are ambiguous (8 if i? > 2.5 is adopted). These MMT 
spectra are shown in Figures [BHHl and their properties 
are listed in Table [TJ 

The spectrum of a i?c~ 22 z ~ 1.6 LEG detected the 
Fe II and Mg II absorption lines, which indicates that 
MMT is sensitive enough to detect luminous LBGs. In 
fact, since the surface density of bright LBGs is low, slit- 
mask instruments are not ideal for the bright end. How- 
ever, the entire SDF can be observed with Hectospec, 
so all '^150 Vauto < 23.0 objects can be simultaneously 
observed. 

2.4. Additional Spectra with Subaru/MOIRCS 

The BzK technique, which identifies galaxies with a 
wide range (old and young, dusty and unreddened) of 
properties, could include objects that would also be clas- 
sified as iV[/l/-dropouts. As a check, cross-matching of 
spectroscopically identified star-forming BzICs with the 
GALEX-SDF photometric catalog was performed. Spec- 
tra of BzKs were obtained on 2007 May 3—4 with Subaru 
using the Multi- Object Infrared Cam era and Spectro- 
graph (MOIRCS: fichikawa et"ani2006[ ). 44 sources were 
targeted and 15 were identified by the presence of Ha 
and [N II] or [O II], [O III] and H/3 emission. One of the 
15 was not in the B-band catalog. Among the 14 ob- 
jects, 7 are also classified as A[/F-dropouts and were not 
previously identified (i.e., LRIS or Hectospec targets). 
This included 5 galaxies at z > 1.5 and 2atz = l — 1.5. 



Their properties are included in Table [TJ Among the 7 
BzKs that did not meet the A^t/V^-dropout criteria, 2 are 
below z = 1.5 and the other five are at high-z. For two 
of the high-z BzKs, one was below the NUV — B = 1.75 
cut because it is faint {V > 25.3), thus not considered 
a Af/F-dropout, and the other missed the B — V = 0.5 
selection by having B — V = 0.53.^^ The other three 
sources have low-z neighboring sources that are detected 
in the NUV, which influences the NUV photometry to 
be brighter. The cause of confusion is due to the poor 
resolution of GALEX, which is discussed further in 14.31 
The detail s of these observatio ns and their results are 
deferred to lHavashi et al.l (|2008l ). 

2.5. Summary of Observations 

In order to probe 1.5 < z < 3 with the Lyman break 
technique, deep (>100 ks) GALEX/ NUV imaging was 
obtained. Spectroscopic observations from Keck and 
MMT independently confirm that most AC/V-dropouts 
(with their UV continuum detected spectroscopically) 
are found to be at 1.5 < z < 2.7. 

A summary of the number of LBGs, stars, and low-z 
interlopers identified spectroscopically is provided in Ta- 
ble [21 Among the spectra targeting A[/V-dropouts (i.e., 
excluding MOIRCS spectra), 53% (30/57) were identi- 
fied, and among those, 63% are at z > 1.5. Including 
seven objects with i? = 2.5 — 3.0, the percentages are 
65% and 62%, respectively. These statistics are improved 
with the final selection criteria discussed in § 13.21 

3. PHOTOMETRIC SELECTION OF Af[/y-DROPOUTS 

This section describes the NUV a.nd optical photomet- 
ric catalogs (§ 13.11) and the methods for merging the two 
catalogs. Then in § 13.21 ^8000 At/V-dropouts are empir- 
ically identified with the spectroscopic sample to refine 
the selection criteria. 

3.1. Revised NUV Photometric Catalogs 

Prior to any measurements, an offset (Aa = —0.39", 
AS = —0.18") in the NUVimage coordinates was applied 
to improve the astrometry for alignment with Suprime- 
Cam data. The scatter in the astrometric corrections 
was found to be (7aq=0.39" and (ta5=0.33". This only 
results in a 0.01 mag correction for A[/F measurements, 
and is therefore neglected. 

The coordinates of ~100000 SDF B-band sources 
with i?auto < 27.5 were used to measure NUV fluxes 
within a 3.39" (2.26 pixels) radius aperture with the 
iraf/daophot task, phot. For objects with iVf/V^ pho- 
tometry below the 3a background limit, the ia value is 
used. This limit is determined from the mode in an an- 
nulus with inner and outer radii of 22.5" and 37.5" (i.e., 
an area of 1200 pixels), respectively. For sources detected 
in the NUV, a point-source aperture correction of a fac- 
tor of «1.83 is applied to obtain the "total" At/I/ flux. 
This correction was determined from the point spread 
function (PSF) of 21 isolated sources distributed across 
the image. The NUV catalog is then merged with the B- 
band catalog from SExtractor (SE; iBertin fc Arnout"i 
Il996f ) that contains BVRci' z' photometry. 

^■^ If the selection criteria were modified to include this object, 
no low-2 interlopers or stars would have contaminated the criteria. 
However, a. B — V < 0.5 is still adopted for simplicity. 
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Fig. 3. — LRIS spectra of confirmed Af?71^-dropouts from the ver. 1 catalog with known redshifts. Most of the LBGs with Lya emission 
are shown here with the remaining in Figure |4] Overlayed on these spectra is the composite template (shown as grey) with the highest 
i?-value (see Table[TJ from cross-correlation. Note that these overlayed templates are intended to show the location of spectral features, and 
is not meant to compare the flux and/or the spectral index differences between the spectra and the templates. The ID number, redshifts, 
and R-values are shown in the upper left-hand corner of each panel. [See the electronic edition of the Journal for a color version of this 
figure.] 



Throughout this paper, "total" magnitudes from the 
Suprime-Cam images are given by SE mag_auto, since 
the corrections between B-band Kron and the 5" diame- 
ter magnitudes were no greater than 0.03 mag for isolated 
(5" radius), point-hke (SE class_star > 0.8) targets. 

The merged catalog was also corrected for galactic ex- 
tinction based on the ICardelli et al.l 1)19891 ) extinction 
law. For the SDF, they are: A{NUV} = 0.137, A{B) 
= 0.067, A{V) = 0.052, A{Rc) = 0.043, A{i') = 0.033, 
and A{z') = 0.025. Since the Galactic extinction for the 
SDF is low, the amount of variation in A{NUV) is no 
more than 0.02, so all A^[/F magnitudes are corrected by 
the same value. 

3.2. Broad-band Color Selection 

Using the sample of spectroscopically confirmed z > 
1.5 LBGs, low-z interlopers, and stars, the color selection 
is optimized to minimize the number of interlopers while 
maximizing the number of confirmed LBGs. In Figure [9l 



known LBGs are identified in the NUV — B versus B — 
V diagram, where the NUV — B color is given by the 
"total" magnitude and the i? — y is the color within a 
2" aperture. The latter was chosen because of the higher 
S/N compared to larger apertures. The final empirical 
selection criteria for the LBG sample are: 

NUV -B> 1.75, (1) 

B-V < 0.50, and (2) 

NUV - B >2.4:{B -V) + 1.15, (3) 

which yielded 7964 A^C/V^-dropouts with 21.90 < V < 
25.30. Among the Hectospec and LRIS spectra, these 
selection criteria included all spectroscopic LBGs and ex- 
cluded 4/5 stars and 4/6 (4/9 with R > 2.5) interlopers. 
Therefore, the fraction of TVf/F-dropouts that are con- 
firmed to be LBGs with the new selection criteria is 86% 
(the R = 2.5 cut implies 79%). Note that while the 
i?-band catalog was used (since the B filter is closer in 
wavelength to the NUV), the final magnitude selection 
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Fig. 4. — Same as Figure [3] but some spectra do not have Lya emission. The strong hne seen in the spectrum of 96927 at ~5570A is 
a sky subtraction artifact, and cosmic rays are seen in the spectra of 94093 (at 3780A), 186254 (at 3325A), and 92942 (at 3990A). These 
features are removed in the cross-correlation process. [See the electronic edition of the Journal for a color version of this figure.] 



was in V, to compare with the rest-frame wavelength 
(«1700A) of ~ 3 LBGs in the i?-band. 

To summarize, a A^f/F-optical catalog was created, and 
it was combined with spectroscopic redshifts to select 
7964 7V[/V-dropouts with NUV - B > 1.75, B - V < 
0.50, NUV - B > 2A{B -V) + 1.15, and 21.90 <V < 
25.30. The spectroscopic sample indicates that 14% of 
A'^t/F-dropouts are definite z < 1.5 interlopers. 

4. CONTAMINATION AND COMPLETENESS ESTIMATES 

Prior to constructing a normalized luminosity function, 
contaminating sources that are not LBGs must be re- 
moved statistically. Section H?T] discusses how foreground 
stars are identified and removed, which was found to be 
a 4—11% correction. Section 321 describes the method 
for estimating low-z contamination, and this yielded a 
correction of 34% ± 17%. These reductions are applied 
to the number of A^t/T^dropouts to obtain the surface 
density of z ~ 2 LBGs. Monte Carlo (MC) realizations 
of the data, to estimate the completeness and the effec- 
tive volume of the survey, are described in § 14.31 The 



latter reveals that the survey samples z « 1.8 — 2.8. 

4.1. Removal of Foreground Stars 

The iGunn fc Strvkeil (|1983[ ) stellar track passes above 
the C/F-dropout selection criteria box (as shown in Fig- 
ure [5]). This poses a problem, as objects that are unde- 
tected in the NUV ca.n be faint foreground stars. A sim- 
ple cut to eliminate bright objects is not sufficient, be- 
cause faint halo stars exist in the SDF (as shown later). 
To reduce stellar contamination, additional photomet- 
ric information from the SExtractor BVRci'z' catalogs 
is used. The approach of creating a "clean" s ample of 
point-like sources, as performed bv iRichmondl (j2005[ ). is 
followed. He used the class_star parameter and the 
difference (6) between the 2" and 3" aperture magni- 
tudes for each optical image. A '1' is assigned when 
the class_star value is 0.90 - 1.00 or 0.10 < S < 0.18, 
and '0' otherwise for each filter. The highest score is 10 
[(1+1) x5], which 2623 V^auto = 21.9 - 26.0 objects sat- 
isfied, and is referred to as "perfect" point-like or "rank 
10" object. These rank 10 objects will be used to define 
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Fig. 5. — Same as Figures |3] and |4] but this shows the z < 1.5 interlopers and galactic stars. [See the electronic edition of the Journal 
for a color version of this figure.] 



the stellar locus, since contamination from galaxies is less 
of a problem for the most point-like sample. Then ob- 
jects with lower ranks that fall close to the stellar locus 
will also be considered as stars after the locus has been 
defined. 

Unfortunately, distant galaxies can also appear 
point-like, and must be distinguished from stars. This 
is done by comparing their broad-band optical colors rel- 
ative to the stellar locus. Figure fTOl shows the B — V, 
V — Rc, and Rc — z' colors used in lRichm ond (2005') for 
the "clean" sample. The stellar locus is defined by the 
solid black lines using brighter (V < 23.0) sources. Fig- 
ure [TO] shows differences in the colors between the stellar 
locus defined f or point-like SDF s tars and those of Gunn- 
Stryker stars. iRichmondl (|2005f) states that this is due 
to metallicity, as the SDF and Gunn-Stryker stars are 
selected from the halo and the disk of the Galaxy, re- 
spectively. 

For each object in the clean sample, the V — Rq color 
is used to predict the B — V and Rq — z' colors along 
the stellar locus (denoted by 'S.L.' in the subscript of the 
colors below). These values are then compared to the ob- 
served colors to determine the magnitude deviation from 
the stellar locus, A = -[(B - T^)obs - {B - F)s.l.] + 
[(i?c — -^Oobs — (-Rc — -^Os.L.]- Therefore, an object with 
A « mag is classi fied as a sta r . Thi s method is similar 
to what is done in 'Richmond (2005) , where an object 
is considered a star if it is located within the stellar lo- 
cus "tube" in multi-color space. This approach provides 
stellar contamination at faint magnitudes , which is diffi- 
cult spectroscopically (jSteidel et al.ll2003[ ). A histogram 
showing the distribution of A in Fig. [TTk reveals two 
peaks: at A « and 0.8 mag. The comparison of A 
versus the y-band magnitude is shown in Fig. Illb . and 
a source is identified as a star if it falls within the se- 
lection criteria shown by the solid lines in this figure. 
A total of 1431 stars V < 26.0 are identified, while the 
remaining 1192 sources are classified as galaxies. The 
surface density as a function of magnitude for the iden- 



tified stars agrees with predictions made bv I Robin et al.l 
(|2003f ) and other surface density measurements near the 
galactic pole. When the A^C/F-dropout selection criteria 
are applied^'^, these numbers are reduced to 336 stars 
(i.e., a 4% contamination for the A^f/l^-dropout sample) 
and 230 galaxies with 21.9 < Vauto < 25.3. 

Sources that are ranked 7 — 9 are also considered and 
were classified as a star or a galaxy using the above ap- 
proach. Of those that met the ATf/F-dropout criteria, 535 
and 252 have the colors of stars and galaxies, respectively. 
Thus, the photometric sample of f/l^-dropouts contains 
7093 objects after statistically removing 871 stars (11% 
of the A^f/l^-dropout) that are ranked 7—10. The reasons 
for only considering objects with a rank of 7 or greater 
are (1) the stellar contamination does not significantly 
increase by including rank 6 or rank 5 objects (i.e., an- 
other 128 rank 6 stars or 1.5% and 143 rank 5 stars or 
1.8%), and (2) comparison of the surface density of rank 
7—10 stars with expectations from models showed ev- 
idence for possible contamination from galaxies at the 
faint end {V > 24.0; A. Robin, priv. comm.), and the 
problem will worsen with rank 5 and 6 objects included. 
As it will be apparent later in this paper, stellar contam- 
ination is small and not expected to significantly alter 
any discussion of differences seen in the luminosity func- 
tion. A hard upper limit by considering objects of rank 
1 and above as stars would imply an additional (rank 1 
to 6) stellar contamination of 14.5%. 

Among the 5 sources spectroscopically determined 
to be stars, 3 of them (71239, 66611, and 149720) are 
classified as stars with the A method, and the other two 
stars (86900 and 178741) fall outside the A selection cri- 
teria. Among the known LBGs, 8 are rank 8—10 and 3 
(166380, 78625, and 133660) are classified as not being 
stars. Since the spectroscopic sample of rank 10 objects 
is small, additional spectra will be required to further 

The B — V and NUV — B color cuts limit the stellar sample 
to spectral types between AO and G8. 
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Fig. 6. — Same as Figures |3] and |4] but these are Hectospec observations of LBGs in the final photometric catalog. The cross-correlation 
template and the typical sky spectrum are shown above and below the spectrum of the source, respectively. [See the electronic edition of 
the Journal for a color version of this figure.] 



optimize the A technique. However, the spectroscopic 
sample (presented in this paper) indicates that 3 — 7% 
of A^f/F-dropouts are stars, which is consistent with the 
4 — 11% derived with the A method. 

4.2. Contamination from z < 1.5 Interlopers 

One of the biggest concerns in any survey targeting 
a particular redshift range is contamination from other 
redshifts. The spectroscopic sample of A^[/F-dropouts 
shows that 5% are definite z < 1.5 galaxies. This num- 
ber increases to an upper value of 51% if the ambiguous 
#[/F-dropouts (that meet the color selection criteria) are 
all assumed to be low-z interlopers. However, it is un- 
likely that all unidentified A'^t/F-dropouts are low-z, since 
LBGs without Lya emission in their spectra^** are likely 
missed. A secondary independent approach for estimat- 
ing low-z contamination, which is adopted later in this 

Either because they do not possess Lya in emission or they 
are at too low of a redshift for Lya to be observed. 



paper, is by using a sample of z < 1.5 emission-line galax- 
ies identified with narrow-band (NB) filters. Si nce a de- 
tailed description of this sample is provided in iLv et all 
((2007 ), only a summary is given below: 
A total of 5260 NB emitters are identified from their ex- 
cess fluxes in the NB704, NB711, NB816, or NB921 filter 
either due to Ha, [O III], or [O II] emission line in 12 
redshift windows (some overlapping) at 0.07 ^ z < 1.47. 
These galaxies have emission line equivalent widths and 
fluxes as small as 20A (observed) and a few xlO~^* erg 
s^^ cm~^, and are as faint as V = 25.5 — 26.0. Cross- 
matching was performed with the A'^t/F-dropout sam- 
ple, which yielded 487 NB emitters as A^f/l^-dropouts. 
The redshift and l^-band magnitude distributions are 
shown in Figure [121 Note that most of the contami- 
nating sources are at 1.0 < z < 1.5, consistent with the 
spectroscopic sample. 

Since this sample represents a fraction of the 0.07 < 
z < 1.5 redshift range, the above results must be inter- 
polated for redshifts in between the NB redshifts. It is 
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Fig. 7. — Same as Figure[6] [See the electronic edition of the Journal for a color version of this figure.] 



assumed that emission-line galaxies exist at all redshifts, 
and possess similar properties and number densities to 
the NB emitters. One caveat of this approach is that 
blue galaxies that do not possess nebular emission lines, 
may meet the A^f/V^-dropout selection. The statistics 
of such objects are not well known, since spectroscopic 
surveys are biased toward emission line galaxies, due to 
ease of identification. Therefore, these contamination es- 
timates are treated as lower limits. A further discussion 
of this approach is provided in § [T] 

Using the redshift distribution shown in Figure [T^ the 
number of objects per comoving volume (N/ AV) is com- 
puted at each NB redshift window. For redshifts not 
included by the NB filters, a linear interpolation is as- 
sumed. Integrating over the volume from z = 0.08 to 
z — 1.5 yields the total number of interlopers to be 
2490 ± 1260, which corresponds to a contamination frac- 
tion of /contam= 0.34±0.17. The error on /contain is from 
Poissonian statistics for each redshift bin, and are added 

Red galaxies are excluded by the B ~V < 0.5 eriterion. 



in quadrature during the interpolation step for other red- 
shifts. This is also determined as a function of magni- 
tude (hereafter the "mag. -dep." correction), since the 
redshift distribution will differ between the bright and 
faint ends. The /contam (V^-band magnitude range) val- 
ues are 0.39 ± 0.20 (22.9-23.3), 0.40 ± 0.21 (23.3-23.7), 
0.37 ± 0.21 (23.7-24.1), 0.31 ± 0.16 (24.1-24.5), 0.27 ± 
0.14 (24.5-24.9), and 0.39 ± 0.19 (24.9-25.3). 

4.3. Modelling Completeness and Effective Volume 

In order to obtain an accurate LF for TVf/V^-dropouts, 
the completeness of the sample must be quantified. 
This is accomplished with MC simulations to calculate 
P{m,z), which is the probability that a galaxy of ap- 
parent y-band magnitude m and at redshift z will be 
detected in the image, and will meet the 7Vf/F-dropout 
color selection criteria. The effective comoving volume 
per solid area is then given by 
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Fig. 9.— NUV - B and S - y colors for 22.0 < Vauto < 25.3 
sources. A total of ~33,000 sources are represented here, but only 
one-third are plotted, for clarity. Sources undetected (at the 3a 
level) in the NUV are shown as grey unfilled triangles while the 
detected sources are indicated as dark grey unfilled squares. Filled 
(unfilled) circles correspond to sources that have been confirmed as 
LBGs with (without) emission lines. Low-z interlopers are shown 
as filled squares while stars are shown as unfilled stars. Skeletal 
stars represent Gunn-Stryker stars. [See the electronic edition of 
the Journal for a color version of this figure.] 



where dV/dz/fl is the differential comoving volume per 
dz per solid area at redshift z. Dividing the number 
of #[/V^-dropouts for each apparent magnitude bin by 
VoS will yield the LF. This approach accounts for color 
selection biases, limitatio ns (e.g., the depth and spatial 
resolution) of the images (jSteidel et al.lll999l ). and choice 
of apertures for "total" magnitude. 

In order to determine P{m,z), a spectral syn- 
thesis model was first constructed from GALAXEV 
(|Bruzual fc Chariot! l2003f) by assuming a constant SFR 
with a Salpeter initial mass function (IMF), solar metal- 
licity, an age of IGyr, and a redshift between z = 1.0 
and z = 3.8 with Az = 0.1 increments. The model 
was reddened by assuming an extinction law following 
iCalzetti et aP feOOO.) with E{B - V) = 0.0 - 0.4 (0.1 in- 
crements) an d modified by accounting for IGM absorp- 
tion following ' Madaul (I1995D. The latt er was chosen over 
other IGM models (e.g.. lBershadv e~a l. 1999) for consis- 
tency with previou s LBG studie s . Th is model is nearly 
identical to that of iSteidel et al.l (|1999i ). 

Figure fT3l shows the redshift evolution of the NUV — B 
and B — V colors for this model. These models were 
scaled to apparent magnitudes of V" = 22.0—25.5 in incre- 
ments of 0.25. These 2175 (29 x 15 x 5) artificial galaxies 
are randomly distributed across the NUV, B, and V im- 
ages with the appropriate spatial resolution (assumed to 
be point-like) and noise contribution with the IRAF tasks 
mkobject (for optical images) and addstar (using the 
empirical NUV PSF). Because of the poor spatial resolu- 
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Fig. 10. — Two color-color diagrams for rank 10 point-like objects. Grey (small) and black (large) squares represent sources brighter than 
Vauto = 26.0 and 23.0, respectively. The Gunn-Stryker stars are shown as stars, and the SDF stellar locus of Richmond (2005) is shown 
as filled squares. The solid lines define the stellar locus for calculating A (see ii 14.11 1. The five sources that have been spectroscopically 
identified to be stars are shown as filled green circles. [See the electronic edition of the Journal for a color version of this figure.] 
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Here fx is the fiux density in the X filter. 

tion of GALEX, each iteration of 435 sources (for a given 
E[B — V] value) was divided into three sub- iterations to 
avoid source confusion among the mock galaxies. The 
artificial galaxies were then detected in the same man- 
ner as real sources. This process was repeated 100 times. 
Note that 21% of artificial sources did not meet the NUV- 
dropout criteria (see e.g., Figure [T4|). as they were con- 
fused with one or more nearby sources detected in the 
NUV. This serves as an estimate for incompleteness due 



to confusion, and is accounted for in the final LF. These 
results are consistent with MOIRCS spectra that finds 
that 14 — 29% of BzKs with 2: > 1.5 was missed by 
A^[/F-dropout selection criteria with nearby objects af- 
fecting the NUV flux. In addition, this simulation also 
revealed that among all mock LBGs with z < 1.5, 30% 
were photometrically scattered into the selection criteria 
of iV[/F-dropouts, which is consistent with the 34% low- 
z contamination fraction predicted in § 14.21 
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Fig. 12. — Redshift (top) and V-band magnitude (bottom) distri- 
butions of 487 NB emitters that meet the Aff7V-dropout criteria. 
Note that the redshift bins are made larger to clearly show the 
histogram. 
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Fig. 13.— Modelled NVV - B and B - V colors for NUV- 
dropouts. The solid, dotted, short-dashed, long-dashed, and dot 
short-dashe d lin es correspond to the spectral synthesis model de- 
scribed in §|43]with E{B -V) = 0.0, 0.1, 0.2, 0.3, and 0.4, respec- 
tively. The thick solid black lines represent the selection criteria in 
§132] 

Figure [14] shows P{m, z) as a function of magnitude 
for E{B -V)= 0.1, 0.2, and 0.0 - 0.4. The latter is de- 
termined from a weighted average w here the E{B — V) 
distribution from'Steidel et al] ()1999D is used for weight- 
ing each completeness distribution. This corresponds to 
an average E{B — V) ~ 0.15. The adopted comoving vol- 
ume uses the weighted-average results. Table [3] provides 
the effective comoving volume per arcmin^, the average 
redshift, the FWHM and standard deviation of the red- 
shift distribution for subsets of apparent magnitudes. 

4.4. Summary of Survey Completeness and 
Contamination 

Using optical photometry, 871 foreground stars (i.e., 
a 11% correction) were identified and excluded to yield 
7093 candidate LBGs. Then z < 1.5 star- forming galax- 



ies, identified with NB filters, were cross-matched with 
the A^f/l^-dropout sample to determine the contamina- 
tion fraction of galaxies at z < 1.5. Redshifts missed by 
the NB filters were accounted for by interpolating the 
number density between NB redshifts, and this yielded 
2490 ± 1260 interlopers, or a contamination fraction of 
0.34 ±0.17. 

To determine the survey completeness, the VcS was 
simulated. This consisted of generating spectral synthe- 
sis models of star-forming galaxies, and then adding ar- 
tificial sources with modelled broad-band colors to the 
images. Objects were then detected and selected as 
A^C/l^-dropouts in the same manner as the final photo- 
metric catalog. These MC simulations predict that the 
survey selects galaxies at z 2.28 ± 0.33 (FWHM of 
z — 1.8 — 2.8), and has a maximum comoving volume of 
2.8 X 10^ Mpc^ arcmin-2_ 

5. RESULTS 

This section provides the key measurements for this 
survey: a z ~ 2 rest-frame UV luminosity function for 
LBGs (§ 15. ip . and by integrating this luminosity func- 
tion, the luminosity and SFR densities are determined 
(§[53). 

5.1. The 1700 A UV Luminosity Function 

To construct a luminosity function, a conversion from 
apparent to absolute magnitude is needed. The distance 
modulus is miroo — Mnoo ~ 45.0, where it is assumed 
that all the sources are at z w 2.28 and the K-correction 
term has been neglected, since it is no more than 0.08 
mag. The luminosity function is given by 



$(M: 



1700 J 



1 ^raw ( 1 /contain ) 

V;ff(Mi7oo) ' 



(5) 



where A^raw is the raw number of A'^f/F-dropouts within 
a magnitude bin (Am = 0.2), K!ff(Afi7oo) is the effective 
comoving volume described in § 14. 3i and /contain is the 
fraction of A^f/l^-dropouts that are at z < 1.5 (see § I4.2|) . 
The photometric LF is shown in Figure[T5] For the mag.- 
dep. /contain case, the adopted correction factor for V < 
22.9 is /contam= 0.34 (the average over all magnitudes). 

Converting the 'Schechter (1976) formula into absolute 
magnitude, the LF is fitted with the form: 
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(6) 

where a; = 10 0-4(a/i7oo J^^itoo). In order to obtain the 
best fit, a MC simulation was performed to consider the 
full range of scatter in the LF. Each datapoint was per- 
turbed randomly 5x10^ times following a Gaussian dis- 
tribution with 1(7 given by the uncertainties in Each 
iteration is then fitted to obtain the Schechter param- 
eters. This yielded for the mag. -dep. /contam case: 
M*7oo = -20.50 ± 0.79, log0* = -2.25 ± 0.46, and 
a = — 1.05±1.11 as the best fit with Icr correlated eriors. 
Since these Schechter parameters are based on lower lim- 
its of low-z contamination (see § 14. 2p . they imply an 
upper limit on (p* . This luminosity function is plotted 
onto Figure [Is] as the solid black line, and the confidence 
contours are shown in F igure [TBI With th e faint-end 
slope fixed to a = -1-60 (jSteidel et al.lll999[ ) and -1.84 
(jReddv et al.ll2008f ). the MC simulations yielded {M^^ooi 
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Fig. 14. — Monte Carlo completeness estimates as a function of redshift for different apparent m agnitude. From left to right is the result 
for E{B — V) = 0.1, 0.2, and 0.0 — 0.4 (a weighted average assuming the E{B — V) distribution of lSteidel et al.lll999l ). [See the electronic 
edition of the Journal for a color version of this figure.] 




-B- This study z = 2,35±0.29 mag, -dep. f^^„, 
-■- Steidel et al- (1999) z=3.04±0.24 
-O- Steidel et al. (1999) z=4.13±0.26 
-•- Reddy et al. (2008) z = 2.20±0.32 BXs 
■A- Rawieki and Thompson (2006) z = 2.20±0.32 BXs 
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Fig. 15. — The observed V-band luminosity function for NUV- 
dropouts. The LF of this work is shown by the thick black 
solid curve with u nfilled squares . Grey points are those excluded 
from the MC fit. ISteidel et al.l II1999I ) measurements are shown 
as filled squares with solid thin curve (z ~ 3) and opened cir - 
cles with short-dashed thin curve {z ~ 4). IReddv et all ll2008l) 
BX results are shown as fille d circles with long-dashed line, and 
ISawicki fc ThompsonI Il2006al) is represented by unfilled triangles 
and dotted line. Corrections to a common cosmology were made 
for Steidel ct ah (1999) measurements, and SFR conversion follows 
IKennicutH | |199^ . [See the electronic edition of the Journal for a 
color version of this figure.] 

log (/)*) of (-20.95±0.29, -2.50±0.17) and (-21.30±0.35, 
-2.75 ± 0.21), respectively. 

5.2. The Luminosity and Star- Formation Rate 
Densities 

The LF is integrated down to Afiyoo — —20.11 — the 
magnitude where incompleteness is a problem — to obtain 
a comoving observed specific luminosity density (LD) of 
log Aim = 26.28 ± 0.69 erg s^^ Hz^^ Mpc'^ at 1700A. 
The conversion between the SFR and specific luminosity 
for 1500-2800Ais SFRuv(Mq yr'^) = 1.4 x 10-28i^(erg 
s^^ Hz~^), where a Salp eter IMF with masses from 
0.1 - IOOMq is assumed ()Kennicuttl [l998[ ). Therefore, 
the extinction- (adopted E[B — V]= 0.15 and Calzetti 
law) and completeness-corrected SFR density of 2; 2 
LBGs is logpstar = -0-99± 0.69 Mq yr"! Mpc'^. Us- 
ing the iMadau et al.l (|1998f) conversion would decrease 






-2.8 -2.6 -2.4 -2.2 

log(0») 

Fig. 16. — Confidence contours representing the best-fitting 
Schechter parameters for the LF. (Top) The mag. -dep. correc- 
tion where the faint-end slope is a free parameter. The vertical 
axes show a while the horizontal axes show log (</>*) (left) and 
M* (right). (Bottom) M* vs. log(0*) for a = -1.6 (left) and 
a = —1.84 (right). The inner and outer contours represent 68% 
and 95% confidence levels. 

the SFR by ~10%. Integrating io L = O. lLt--^, where 
is at z ~ 3 (A4^^3 = -21.07, ISteidel et all 
I1999I ). yields log£ = 26.52 ± 0.68 erg s'^ Hz"! Mpc'^ 
or an extinction-corrected SFR density of log pstar = 
-0.75 ± 0.68 Mq yr-i Mpc'^.i^ 

5.3. Summary of Results 

A UV luminosity function was constructed and yielded 
a best Schechter fit of A/j^^oo = -20.50 ± 0.79, log(/-* = 
-2.25±0.46, and a = -1.05±1.11 for z - 2 LBGs. The 
UV specific luminosity density, above the survey limit, is 
log Aim = 26.28 ±0.68 erg s^^ Hz"! Mpc'^. Correcting 
for dust extinction, this corresponds to a SFR density of 
logp.tar = -0.99 ± 0.68 Mq yr"! Mpc'^. 

6. COMPARISONS WITH OTHER STUDIES 

Comparisons in the UV specific luminosity densities, 
LFs, and Schechter parameters can be made with pre- 
vious studies. First, a comparison is made between the 

The above numbers are upper limits if the lo w-z c ontamina- 
tion fraction is higher than estimates described in § 14.21 
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z ~ 2 LBG LF with z - 2 EX and z - 3 LBG LFs. 
Then a discussion of the redshift evolution in the UV lu- 
minosity density and LF (parameterized in the Schechter 
form) is given in § 16.21 

The results are summarized in Figures [151 UHl ^^nd [19] 
and Table H) For completeness, three different UV spe- 
cific luminosity densities are reported by integrating the 
LF down to: (1) 0.1L*^3; (2) Lum, the limiting depth of 
the survey; and (3) L — 0. The latter is the least confi- 
dent, as it requires extrapolating the LF to the faint-end, 
where in most studies, it is not well determined. 

6.1. UV-selected Studies at z ~ 2 — 3 

In Figure [l5l the z - 2 LBG L F at the bright end 
is similar to those of LBGs from Stcidcl ct al. (19991) 
and BX galaxies fro m iSawicki fc Thompson (20()6a) and 
iReddy et al.l ()2008f ): however, the faint end is systemat- 
ically higher. This is illustrated in Figure [T7I where the 
ratios between the binned z ^ 2 UV LF a nd the fitted 
Schech ter forms of lSteidel et al. (1999) and R eddv et al.l 
()2008[ ) are shown. When excluding the four brightest 
and two faintest bins, the A^f/V^-dropout LF is a factor 
of 1-7 =b 0.1 with respect to z ~ 3 L BGs of ISteidel et al.l 
(119991) and z ^ 2 BX ga laxies _ of i Reddv et al.l (|2008l ) 



and lSawicki fc Thompson! (|2006a|) . The hard upper limit 
for stellar contamination (see § 14. ip would reduce this 
discrepancy to a factor of 1 .4 ± 0.1. There app ears to 
be a trend that the ratio to iReddv et all (|2008f) LF in- 
creases towards brighter magnitudes. This is caused by 
the differences in the shape of the two LFs, particu- 
larly the faint-end slope. The in crease in the r atio i s 
less noticeable when compared to ISteidel et al.l (|1999D . 
w hich has a shallower f aint-en d slo pe. Since the LFs 
of lSawicki fc Thompson! ^006^ and !Reddv eFal! (!2008! ) 
are similar, the comparison of any results between the 
A^C/F-dropout and the BX selections will be made di- 
rectly against !Reddv et all (!2008! ). 

All 11 points are 1 — 3a from a ratio of 1. It has been 
assumed in this comparison that the amount of dust ex- 
tinction does not evolve from z ~ 3 to z ~ 2. Evidence 
supporting this assumption is: in order for the intrinsic 
LBG LFs at z 2 and 3 to be consistent, the popula- 
tion of LBGs at z ~ 2 would have to be relatively less 
reddened by AE{B - V) ^ 0.06 (i.e., E[B - V] ^ 0.09 
assuming a Calzetti extinction law). However, the stel- 
lar synthesis models, described previously, indicate that 
E{B — V) ^ 0.1 star-forming galaxies are expected to 
have observed B - V ^ 0.1, and only 15% of NUV- 
dropouts have B ~ V < 0.1. This result implies that 
dust evolution is unlikely to be the cause for the discrep- 
ancy seen in the LFs. 

To compare the luminosity densities, the binned LF 
is summed. This is superior to integrating the Schechter 
form of the LF as (1) no assumptions are made between 
individual LF values and for the faint-end, and (2) the 
results do not suffer from the problem that Schechter 
parameters are affected by small fiuctuations at the 
bright- and faint-ends. The logarithm of the binned lu- 
minosity densities for —22.91 < Mimn < —20.11 ar e 
26.27±0.16 (this work), 26.02±0.04 (ISteidel et al."1999D. 
and 26.08 ± 0.07 ergs s^^ Rz'^ Mpc'^ ([Reddv et al] 
120081 ) ■ which impfics that the z ~ 2 LBG UV luminos- 
ity density is 0.25 ± 0.16 dex higher than the other two 
studies at the 85% confidence level. 




Fig. 17. — Comparisons of the LBG LF with ot her LFs. The 
rati os of the z ~ 2 LBG L F to the Schechter fits of lSteidel et all 
((1999) LF and R eddv et al.1 1| 2008 ) are shown in the top and bottom 
panels, respectively. On average, the z ~ 2 LBG LF is a factor of 
1.7 ± 0.1 higher than these studies. 

Since the low-z contamination fraction is the largest 
contributor to the errors, more follow-up spectroscopy 
will reduce uncertainties on the LF. This will either con- 
firm or deny with greater statistical significance that the 
luminosity density and LF of z ~ 2 LBGs are higher than 
the z ~ 3 LBGs and z ~ 2 BXs. 

6.2. Evolution in the UV Luminosity Function and 
Density 

The Schechter LF parameters, listed in Table (H are 
plotted as a function of redshift in Figure [TH There 
appears to be a systematic trend that M* is less neg- 
ative (i.e., a fainter L*) by «1 mag at higher redshifts 
for surveys with a < —1.35. No systematic evolution is 
seen for 0*, given the measurement uncertainties. Lim- 
ited information are available on the faint-end slope, so 
no analysis on its redshift evolution is provided. It is of- 
ten difficult to compare Schechter parameters, since they 
are correlated, and without confidence contours for the 
fits of each study, the apparent evolution could be in- 
significant. A more robust measurement is the product 
{(/)* X L*), which is related to the luminosity density. 

The observed LDs, integrated to 0.1L*^3, show a slight 
increase of « 0.5 dex from z ~ 6 to z ~ 3. How- 
ever, the two other luminosity densities appear to be 
fiat, given the scatter in the measurements of « 0.5 — 1.0 
dex. A comparison between z ~ 2 and z ^ 5 stud- 
ies reveal a factor of 3 — 6 higher luminosity density at 
z ~ 2. The extinction-corrected results for = 

and = O.lLt-q show a. factor of 10 increase from 
6 iBouwens et al.l (!2007D 's measurement to z 



IBouwens et al.l (!2007D assumed a lower dust extinction 
correction. If an average E{B — V) = 0.15 with a Calzetti 
law is adopted, the rise in the extinction-corrected lumi- 
nosity density is « 3. 

7. DISCUSSION 

In this section, the discrepancy between the UV LF of 
this study and two BX studies, shown in § 16.11 is exam- 
ined. Three possible explanations are considered: 



Lyman Break Galaxies at z = 1.8 — 2.8 



15 



-e- 



o 



-3.5 



-21.5 



-21 



I I 1 1 I I 1 1 I 1 1 I I 1 1 I 1 1 I I 1 1 I I III I 1 1 I I I I I 1 1 I I I I I I 1 1 I I I I I III 1 1 III 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



— hE3h 



I I I I I I I I I I I I I I I I I I 



— 



hE3h 



< a t: -!.35 
I I I I I I I I I I I I I I I I I I I I I 



a > -11.35 



22 M I I 1 I I I I I I 1 I I I I I I I I I I I I I Ml I I I I I I I I I I I I I I I I I I I I I I I III 1 I I 1 I 1 I I 1 I I I I 1 I I I I 1 I I I 



-20.5 - 



-20 



I I I I I I I I I I I I I I I I I I 



e 



-1.70 < a ^ -1.35 

I I I I I I I I I I I I I I I I I I I I I 



B 



a > -1.35 



2 3 4 5 
redshift 







■ 


□ O Ly+(2008) 






A 


Iwata + (2007) 


A 


ST(2006) 


A 


Yoshida+(2006) 


• 


Hildebrandt+(2007) 


• 


Reddy+(2008) 


• 


Wadadekar+(2006) 


o 


Giavalisco + (2004) 


o 


Paltani + (2007) 


O 


Steidel + (1999) 


■ 


Burgarella + (2007) 


□ 


Massarotti + (2001) 


■ 


Steidel + (1999) 


□ 


Bouwens + (2006,2007) 




Madau + (1996) 


□ 


Shim+(2007) 



Fig. 18. — Compiled Schechter parameters of LBG and BX studies versus redshift. Top and bottom show the the normalization {</>*), 
and the "knee" of the UV LF {M*), respectively. Measurements are grouped according to a: < —1.70, between —1.70 and —1.35, and 
> —1.35. This A'^f/y-dropout work is show n as black filled square (a = — 1.05). The color and symbol conventions for studies in Figure [Tsl 
are identical for this figure. In the legend. ISawicki fc ThompsonI Il2006al ) is abbreviated as "ST(2006)". Some points are not shown here 
but have luminosity density measurements presented in Figure [Tol 

1. Underestimating low-z contamination. To es- 
timate contamination, a large sample of z < 1.5 NB 
emitters was cross-matched with the A'^f/F-dropout sam- 
ple. This method indicated that 34% ± 17% of NUV- 
dropouts are at z < 1.5. However, it is possible that 
star- forming galaxies at z = 1 — 1.5 could be missed 
by the NB technique, but still be identified as NUV- 
dropouts. This would imply that the contamination rate 
was underest imated. To shif t the iV?7V^-dropout LF to 
agree w ith iRe ddv et al.1 ()2008f) and lSawicki fc ThompsonI 
(|2006al) would require that the contamination fraction 
be more than 60%. However, the spectroscopic sam- 
ple has yielded a large number of genuine LBGs and a 
similar low-z contamination (at least 21% and at most 
38%). If the large (60%) contamination rate is adopted, 
it would imply that only 15 of 40 spectra (LRIS and 
Hectospec) are at z > 1.5, which is argued against at 
the 93% confidence level (98% with R = 2.5 threshold), 
since 24 LBGs (1.6 times as many) have been identified. 
Furthermore, the LRIS and Hectospec observations in- 
dependently yielded similar low contamination fractions, 
and the MC simulation (that involved adding artificial 
LBGs to the images) independently suggested 30% con- 
tamination from z < 1.5. 



2. Underestimating the comoving effective vol- 
ume. The second possibility is that Vcs was underes- 
timated, as the spectral synthesis model may not com- 
pletely represent the galaxies in this sample, and misses 
z ~ 1 — 1.5 galaxies. However, a comparison between a 
top-hat P{m, z) from z = 1.7 — 2.7 versus z = 1.4 — 2.7 
(z = 1.0 — 2.7) would only decrease number densities by 
« 20% (37%). Note that the latter value is consistent 

with /cont am- 

3. Differences between LBG and BX galaxies se- 
lection. This study uses the Lyman break technique 
while other studies used the 'BX' method to identify 
z ~ 2 galaxies. Because of differences in photometric 
selection, it is possible that the galaxy population iden- 
tified by one method does not match the other, but in- 
stead, only a fraction of BX galaxies are also LBGs and 
vice versa. This argument is supported by the higher 
surface density of LBGs compared to BXs over 2.5 mag, 
as shown in Figure [20k. However, their redshift distribu- 
tions, as shown in Figure [20b . are very similar. 

This scenario would imply that there is an increase in 
the LF and number density of LBGs from z '--^ 3 to z '--^ 2, 
indicating that the comoving SFR density peaks at z ~ 2, 
since there is a decline towards z; ~ from UV studies 
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Fig. 19. — The observed (left) and extinction-corrected (right) UV specific luminosity densities as a function of redshift. The luminosity 
function is integrated to three different limits: L = (top panel), L = (the survey's limit; middle panel) , and L = . 1L*_ ci. The color 

and point-type schemes are the same as Figure [TH] The SFR densities are shown on the right axes following iKennicuttI (|l99f) conversion. 
For the 2 ~ 2 LBG luminosity density integrated to L = Liim, only one value is shown, since all the fits with different a are almost identical. 

fication, 871 foreground stars have been identified 
and removed from the photometric sample. This 
corresponds to a 4 — 11% correction to the NUV- 
dropout surface density, which is consistent with 
the 3 — 7% from hmited spectra of stars presented 
in this paper. 

4. In addition, low-z contamination was determined 
using a photometric sample of NB emitters at 
z 5, 1.47. This novel technique indicated that 



(see lHopking|l200^ . and references therein). However, it 
might be possible that the selection (NUV—B — V) of 
z ~ 2 LBGs could include more galaxies than the UnGR 
color selection used to find z ~ 3 LBGs. Although no 
reason exists to believe that z ^ 3 LBG selection is more 
incomplete than at z ~ 2 (nor is there any evidence for 
such systematic incompleteness for z > 4 LBGs), it is 
difficult to rule out this possibility for certain. But if so, 
then the SFR density might not evolve. In addition, the 
conclusion that z 2 is the peak in star-formation is 
based on UV selection techniques, which are less sensi- 
tive at identifying dusty {E[B — V] > 0.4) star-forming 
galaxies. However, spectroscopic surveys have revealed 
that the sub-mra gal axy population peaks at z w 2.2 
(jChapman et al.l[2005f ). which further supports the above 
statement that z 2 is the epoch of peak star-formation. 

8. CONCLUSIONS 

By combining deep GALEX/ NUV and optical 
Suprime-Cam imaging for the Subaru Deep Field, a large 
sample of LBGs at z ~ 2 has been identified as NUV- 
dropouts. This extends the popular Lyman break tech- 
nique into the redshift desert, which was previously dif- 
ficult due to the lack of deep and wide-field UV imaging 
from space. The key results of this paper are: 

1. Follow-up spectroscopy was obtained, and 63% 
of identified galaxies are at z = 1.6 — 2.7. 
This confirms that most iVf/l^-dropouts are LBGs. 
In addition, MMT/Hectospec will complement 
Keck/LRIS by efficiently completing a spectro- 
scopic survey of the bright end of the LF. 

2. Selecting objects with NUV - B> 1.75, B-V < 
0.5, and NUV-B > 2.4(B-\/) + l. 15 yielded 7964 
Aff/y-dropouts with V = 21.9 - 25.3. The spec- 
troscopic sample implied that 50—86% of NUV- 
dropouts are LBGs. 

3. Using broad-band optical colors and stellar classi- 



the contamination fraction is (at least) on average 
34% ± 17%, which is consistent with the spectro- 
scopic samples and predictions from MC simula- 
tions of the survey. 
5. After removing the foreground stars and low-z in- 
terlopers, MC simulations were performed to es- 
timate the effective comoving volume of the sur- 
vey. The UV luminosity function was constructed 
and fitted with a Schechter profile with M^yoq — 
-20.50 ± 0.79, logc/)* = -2.25 ± 0.46, and a = 
-1.05 ± 1.11. 

A compilation of LF and SFR measurements for 
UV-selected galaxies is made, and there appears to 
be an increase in the luminosity density: a factor of 
3—6 (3 — 10) increase from z 5 (z ~ 6) to z ~ 2. 
Comparisons between j V^/F-d ropouts with LBGs 
at z ^ 3 (Steidel et all 1 19991) and BXs at z ~ 2 
dSawicki & Thom psoiJ " l2006a[ iReddv et all l2008f l 
reveal that the LF is 1.7 ±0.1 (1.4 ±0.1 if the hard 
upper limit of stellar contamination is adopted) 
times higher than these studies. The summed lu- 
minosity density for z ~ 2 LBGs is 1.8 times higher 
at 85% confidence (i.e., 0.25 ± 0.16 dex). 
Three explanations were considered for the discrep- 
ancy with z ^ 2 BX studies. The possibility of 
underestimating low-z contamination is unlikely, 
since optical spectroscopy argues against the possi- 
bility of a high (60%) contamination fraction at the 
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Fig. 20. — Surface densities and redshift distributions for 2 ~ 2 BXs and LBGs. In (a), the surface densities of LBGs and BXs are shown 
as circles and triangles, respectively. Both studies have stellar and low-2 contamination corrections applied. This figure reveals that the 
LBG surface density is systematically higher than the BX's. The redshift distributions are shown in (6). T he shaded (unshaded) histogram 
corresponds to BXs (LBGs). For the BX, the redshift distribution is obtained from lReddv et ahl 120081) spectroscopic sample, while the 
LBG is determined from the MC simulations described in § 4.3 for all magnitudes. The similarities in redshifts surveyed by both studies 
and the higher surface density of LBGs indicate that the BX technique misses a fraction of LBGs. 



93% confidence. Second, even extending the red- 
shift range to increase the comoving volume is not 
sufficient to resolve the discrepancy. The final pos- 
sibility, which cannot be ruled out, is that a direct 
comparison between BX-selected galaxies and LBG 
is not valid, since the selection criteria differ. It is 
likely that the BX method may be missing some 
LBGs. This argument is supported by the simi- 
lar redshift distribution of BXs and LBGs, but the 
consistently higher surface density of LBGs over 
2.5 mag. 

9. If the latter holds with future reduction of low-z 
contamination uncertainties via spectroscopy, then 
the SFR density at z ~ 2 is higher than z > 3 and 
z < 1.5 measurements obtained v ia UV selection. 
Comb ined with sub-mm results (iChapman et al.) 
I2005D . it indicates that 2 2 is the epoch where 
galaxy star-formation peaks. 
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APPENDIX 

A. INDIVIDUAL SOURCES OF SPECIAL INTEREST 

In most cases, the confirmed LBGs showed no unique spatial or spectral properties. However, 3 cases are worth 
mentioning in more detail. 

1. SDFJ132431. 8+274214. 3 (179350). Upon careful examination of the 2-D spectra, it appears that the Lya emis- 
sion from this source is offset by wl.l" (9 kpc at 107° east of north) from the continuum emission, which is shown in 
Figure [2Tk. The extended emission appears in the individual exposures of 15 — 30 minutes. The deep (3cr — 28.45) 
i?-band image (Figure [2Tb ) reveals that there are no sources in this direction and at this distance, assuming that the 
continuum emission in the spectrum corresponds to the bright source in the _B-band image. The two sources located 
below the bright object in Figure [2Tb are too faint for their continuum emission to be detected with LRIS. Also, 
absorption features seen in the 1-D spectra (see Figure [3^) are at nearly the same redshift as Lya. This indicates that 
the Lya emission is associated with the targeted s ource, rather than a secondary nearby companion. 

Extended Lya emission galaxies are rare (e.g., [Saito et al.l [20061 . have the largest sa mple of 41 obje cts), and the 
extreme cases are extended on larger 100 kpc) scales, such as LABI and LAB2 of Stei del et al.[ (|200Clf). In addition, 
exten ded Lya emission has been seen in some cases that show evidence for energetic galactic winds ( Mas-Hesse et al.l 
[2003D . Either this source is a fortuitous discovery from a dozen spectra, or perhaps a fraction of A'^t/F-dropouts have 
extended Lya emission. The physical significance of this source is not discussed here, given limited information. 
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Fig. 21. — Optical images for 179350. (a) The 2-D spectrum with wavelength increasing to the right shows hya emission offset by 
Ril" from the center of the continuum. The vertical white line corresponds to 2". (b) The Suprime-Cam _B-band image centered on the 
targeted source shows that there are no sources in the direction of the extended emission. The two white vertical lines correspond to the 
slit, so (b) is rotated to have the same orientation as (a), and the vertical scales are the same. [See the electronic edition of the Journal 
for a color version of this figure.] 




2. SDFJ132452. 9+272128. 5 (62056). The 1- and 2-D spectra for this source reveal an asymmetric emission 
hne, as shown in Figure 122b ,. but with a weak "bump" about lOA blue-ward from the peak of Lya emission. The 
i?-band image (see Figure [53]) shows two nearby sources where one is displaced «2" nearly in the direction of the slit 
orientation while the other source is displaced in the direction perpendicular to the slit orientation. It may be possible 
that the blue excess is originating from the latter source due to a slight misalignment of the slit to fall between the 
two sources (i.e., they are physically near each other). To confirm this hypothesis, spectroscopy with a 90° rotation of 
the slit would show two sources with Lya emission «800 km s~^ apart. 

3. SDFJ132450. 3+272316. 24 (72012). This object is not listed in TablefH as it was serendipitously discovered. 
The slit was originally targeting a narrow-band (NB) emitter. The LRIS-R spectrum showed an emission line at 
7040A, but the blue-side showed a strong emission line that appears asymmetric at w4450A. One possibility is that 
the 4450A feature is Lya, so that the 7040A emission line is the redshifted C III] A1909, but at 2; = 2.6634, C III] is 
expected at ss6994A. This «40A difference is not caused by poor wavelength calibration, as night sky and arc-lamps 
lines are located where they are expected in both the blue and red spectra. In Figure the _B-band image reveals 
two sources, one of which is moderately brighter in the NB704 image, as expected for a NB704 emitter. These two 
sources were too close for SExtractor to deblend, but the coordinate above has been corrected. Because the NB704 
emitter is a foreground source, the measured NUV&ux for the other source is affected, and results in a weak detected 
source in the NUV. Thus, this source is missed by the selection criteria of the ver. 1 catalog and those described in 
§ 13.21 It is excluded from the spectroscopic sample discussed in § [2] 

This source is of further interest because it also shows a blue excess bump (shown in Figure [22]) much like 62056, but 
weaker. This blue bump does not correspond to a different emission line with the same redshift as the 7040A emission 
line. Since the bump is lOA from the strong Lya emission, it is likely associated with the source producing Lya. 
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Fig. 23. — The _B-band image cropped to 20" on a side and centered on 62056. The white box with thick Hnes is the LRIS sht intended 
to target the bright object. However, a 1.5" offset of the slit in the north-west direction (as shown by the thin white box) may explain the 
blue excess seen in the 1-D and 2-D spectra (Figure |22j by including both objects. 

Both 62056 and 72012 were obtained on the second mask. These blue bumps are not due to a misalignment of single 
exposures when stacking the images together, as other equally bright sources in the mask with e mission lines do not 
show a secondary blue peak. Other studies hav e also seen dual peak Lya emission profiles (e.g.. iTapken et ani2004 
l2007t iCooke et al.ll2008t iVerhamme et al]|2008D. In additi on, high resolution spectra of 9 LBGs have also revealed 3 
cases with double-peaked Lya profile ( Shaplev et al.|[2C)06f ). which indicates that such objects may not be rare. 




Fig. 24. — Postage stamp images (10" on a side) for 72012. From left to right is NW, B, Rc, and NB704. North is up and east is to 
the left. The source on the right shows a weak excess in NB704 relative to the broad-band images. 
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TABLE 1 

Properties of Spectroscopically Targeted Af[/V-DROPOUTS and BzKs 
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Name (SDF) 



UV and Optical measurements 



Spectroscopic measurements 
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2 
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2 
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6 
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2 


97 


5 


78733H 


132346.4+272426.09 


4.126 


0.367 


27.164 


23.038 


22.683 


22.586 


22.567 


22.575 


1.8112'' 


2 


15 


7 


66488H 


132520.4+272220.46 


>4.595 


0.447 


>27.165 


22.570 


22.126 


21.960 


21.897 


21.833 


2.4233'' 


2 


89 


5 




1 "^9^1 Q^97441 7 9f\ 


3.652 


0.409 


26.211 


22.559 


22.148 


21.999 


21.969 


21.930 


2.3434'' 


2 


84 


g 






















U.oyoz 


A 


vo 
( y 


[3] 


190947H 


132346.0+274419.92 


>4.344 


0.212 


>27.140 


22.796 


22.597 


22.423 


22.259 


22.081 


1.7273'' 


2 


89 


3 






















0.9612 


2 


38 


[3] 


153628H 


132514.5+273743.52 


>4.563 


0.457 


>27.268 


22.705 


22.243 


22.066 


21.993 


21.962 


0.0551'' 


2 


63 


[2] 










Undetected A f V-dropouts 














174747L 


132436.7+274129.12 


1.951 


0.247 


26.332 


24.381 


24.184 


24.080 


24.028 


23.969 










182447L 


132429.1+274249.80 


>2.603 


0.385 


>27.105 


24.502 


24.126 


23.642 


23.192 


22.584 










180088L 


132421.6+274223.22 


>2.720 


0.129 


>27.120 


24.400 


24.262 


24.198 


24.040 


23.766 










172253L 


132414.3+274100.26 


>2.911 


0.224 


>27.120 


24.209 


24.017 


23.970 


23.960 


23.853 










184387L 


132414.6+274308.58 


>2.362 


0.173 


>27.137 


24.775 


24.606 


24.275 


24.109 


23.586 










63360L 


132433.2+272142.21 


>2.353 


0.443 


>27.167 


24.814 


24.275 


23.840 


23.600 


23.346 










113109L 


132514.6+273028.10 


>2.478 


0.209 


>27.158 


24.680 


24.417 


24.408 


24.252 


24.298 










94367L 


132459.1+272709.00 


>3.386 


0.246 


>27.169 


23.783 


23.567 


23.516 


23.436 


23.309 











Note. — Identified sources are based on an i? > 3.0 criterion (exceptions are AGNs, stars, and those with emission lines). Col. (1) is the S-band catalog ID, Col. (2) is the J2000 coordinates, 
and magnitudes and colors are given in Cols. (3) to (10). The cro ss-correlated redshif ts and fi-values from xcsao are provided in Cols. (11) and (12). The t emplate yieldin g the highest i?-value 
is given in Col. (13), where 1-4 correspond to the four spectra of ISteidel et al 1j200 ^ from stronge st Lya absorption to emission, and 5 and 6 refer to the lShapley et al.l (2003) composite and 
cB5S spectr a, respectively. '7 ' corresponds to the rest-frame NUV spectra from fS avaglio et al.l I I2004D , and 'AGN' refers to a SDSS QSO template. For interlopers, the six SDSS composite spectra 
presented in |Yip et al.| j2004D correspond to [1] to [6] from strongest absorption-line to strongest emission-line systems. [7], [8], and [9] correspond to the RVSAO templates "femtemp", "EA" , and 
"eatemp" , respectively. For objects with Lya emission, the line flux (in units of 10"^^ ergs s~^ cm~^) and rest-frame EWs (in units of A) arc given in Cols. (14) and (15), respectively. These 
were measured using the splot routine. 

a b 

The character following the ID number corresponds to the spectrograph used: 'H'— Hectospcc, 'L'— LRIS, 'M'— MOIRCS. Values in superscript correspond to other templates, which yielded 

similar cross-correlated velocities with R > 2.5.^ As discussed in § ^] this source shows an unusual Lya profile, so the reported flux and rest-frame EW excluded the blue excess by deblending in 

IRAF splot. The values reported here correspond to the best cross-correlated results, but may be wrong due to the low S/N of the spectra. 
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TABLE 2 

Summary of Spectroscopic Observations 



Instrument 


Total 


LBGs [AGNs] 


z < 1.5 


stars 


Ambiguous 


Undetected 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


LRIS 


36 (28) {4} 


12 {2} 


5 (1) {2} 


3(0) 


8 (7) {-4} 


8(8) 


Hectospec 


21 (20) {3} 


7 [2] {2} 


1 (1) {1} 


2 (1) 


11 (11) {-3} 


0(0) 


MOIRCS 


44 


10 (5) 


4(2) 









Total 


101 


29 (24) [2] {4} 


10 (4) {3} 


5(1) 


19 (18) 1-7} 


8(8) 



Note. — Sources with z > 1.5 arc classified as "LBG". Values in square brackets are those that appear to be AGNs, and those in parentheses 
meet the final selection criteria in ^ |3.2| Values in curly brackets represent LBGs that are reclassified as "identified" if a lower (R — 2.5) threshold 
is adopted rather than a. R — 3.0 cut. None of the LBGs and AGNs was missed by the final selection criteria. 



TABLE 3 

Effective Volume Estimates 



V mag 


^avg 




FWHM{z) 








FWHM{z) 




^avg 




FWHM{z) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 






E{B 


~V) = 0.0 






E{B 


-V) = 0.1 








• -V) = 0.2 




22.0-22.5 


2.422 


0.312 


1.860-2.962 


2.93 


2.313 


0.318 


1.757-2.860 


2.89 


2.150 


0.304 


1.627-2.667 


2.69 


22.5-23.0 


2.420 


0.312 


1.865-2.961 


2.67 


2.304 


0.323 


1.746-2.859 


2.67 


2.148 


0.307 


1.613-2.668 


2.48 


23.0-23.5 


2.418 


0.317 


1.862-2.962 


2.36 


2.304 


0.328 


1.738-2.864 


2.40 


2.146 


0.312 


1.592-2.664 


2.23 


23.5-24.0 


2.401 


0.333 


1.831-2.969 


2.12 


2.290 


0.339 


1.717-2.863 


2.20 


2.166 


0.320 


1.585-2.682 


1.97 


24.0-24.5 


2.379 


0.361 


1.772-2.974 


2.14 


2.261 


0.355 


1.645-2.846 


2.10 


2.157 


0.331 


1.570-2.684 


1.93 


24.5-25.0 


2.312 


0.369 


1.680-2.910 


2.14 


2.199 


0.349 


1.575-2.774 


2.02 


2.104 


0.314 


1.561-2.631 


1.81 


25.0-25.5 


2.220 


0.352 


1.603-2.788 


1.80 


2.133 


0.325 


1.569-2.659 


1.44 


2.043 


0.292 


1.560-2.487 


1.01 



-V) = 0.3 

1.563-2.569 

1.573-2.567 

1.573-2.574 

1.567-2.572 

1.559-2.559 

1.557-2.422 

1.550-2.004 



E{B -V) = 0.4 



E{B - y) = 0.0 - 0.4 



22.0-22.5 
22.5-23.0 
23.0-23.5 
23.5-24.0 
24.0-24.5 
24.5-25.0 
25.0-25.5 



2.074 
2.079 
2.088 
2.077 
2.060 
2.002 
1.898 



E[ 
0.290 
0.287 
0.291 
0.295 
0.293 
0.279 
0.266 



B 



2.41 


1.956 


0.236 


1.553- 


-2.359 


1.78 


2.278 


0.333 


1.752- 


-2.837 


2.79 


2.16 


1.952 


0.233 


1.559- 


-2.352 


1.57 


2.274 


0.335 


1.740- 


-2.834 


2.56 


1.88 


1.957 


0.230 


1.562- 


-2.350 


1.40 


2.273 


0.338 


1.741- 


-2.832 


2.28 


1.79 


1.952 


0.232 


1.557- 


-2.349 


1.34 


2.268 


0.345 


1.711- 


-2.839 


2.07 


1.77 


1.934 


0.235 


1.553- 


-2.313 


1.28 


2.248 


0.360 


1.610- 


-2.834 


2.02 


1.29 


1.887 


0.237 


2.131- 


-2.243 


0.56 


2.194 


0.354 


1.574- 


-2.762 


1.91 


0.35 


1.967 


0.334 


1.962- 


-2.062 


0.08 


2.139 


0.336 


1.571- 


-2.652 


1.30 



Note. — Results for MC simulations described in %. |4.3| for different assumed E{B — V) values. Col. (1) lists the apparent V-band magnitude. 
Cols. (2), (6), and (10) show the average redshift {^Zang)-, and Cols. (3), (7), and (11) list the redshift Icr uncertainties. Cols. (4), (8), and (12) 
give the FWHM of the redshift distribution, and the effective comoving volume per area in units of 10^ ^7(f Mpc^ arcmin"^ are in Cols. (5), (9), 
and (13). 



TABLE 4 

Compilation of UV Luminosity Functions and Luminosity Densities 



Reference 


z 


N 


Area 


Arcst 


Ce 


Cl 




log (p* 




a 




log Cobs 






(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


^ U. 11^2 — 3 




L > 




(IS) 


(14) 


Bouwens et aL (2006) 


5.90±0.30 


506 


344.2 


1350 


1.51 


1.000 


1.000 


-2.69«i? 


-20.25±0.20 


-1 70+0.21 
-1. (ci_0 20 


26.047 


26.224l[5;"3 


26.567 


Bouwens et aL (2007") 


3.80±0.35 


4671 


347 


1600 


2.692 


1.000 


1.000 


-2.89tH? 


-20.98±0.10 


-1.73±0.05 


26.276 


26.503±0.05 


26.668 




5.00±0.35 


1416 


367 


1600 


1.995 


1.000 


1.000 


-3.00+°" 

— U. ib 


-20.64±0.13 


-1.66±0.09 


25.952 


26.161±0.06 


26.313 




5.90±0.30 


627 


396 


1350 


1.51 


1.000 


1.000 


9 ocr+0.16 
-^•»O_0.15 


-20.24±0.19 


-1.74±0.16 


25.885 


26.100±0.08 


26.421 


Burearella et aL (2007) 


l.iUltU.ZU 


A on 


0/1 v 

y4 / 


1 cnn 
ioUU 


O.DOl 


1 nnn 


1 nnn 

i.UUU 














Foucaud et aL (2003) 




1294 


1700 


1900 




1.000 


1.000 














Giavalisco et aL (2004) 


3.78±0.34 


1115 


316 


1500 


7.48 


1.000 


1.000 






-1.60f 




26.212±0.07 






4.92±0.33 


275 


316 


1500 


7.48 


1.000 


1.000 






-1.60f 




26.017±0.14 






5.74±0.36 


122 


316 


1500 


7.48 


1.000 


1.000 






-1.60' 




26.061±0.19 




Hildebrandt et aL (2007) 


2.96±0.24 


14283 


9.99 


1650 


3.88 


1.000 


1.000 


-3.29±0.08<= 


-22.43±0.11° 


-1.60'='' 


26.362 


26.097 


26.492 


Jwata et aL (2007) 


4.80±0.40 


853 


1290 


1500 


7.48 


1.000 


1.000 


-3.39«;1 


-21.23±0.30 


-1-^0-0.32 


25.84110:04 


25.679 


25.995lE;:g3 


NIadau et aL ( 1996) 


2.75±0.75 


69 


4.65 


1620 


3.23 


1.108 


0.712 




<26.101 


4 nn+n 50 


14 


4.65 


1630 


3.23 


1.170 


0.662 










^25 588 




Massarotti et aL (2001) 


1.50±0.50 


315 


5.31 


1500 


8.10 


0.980 


0.844 










26.240±0.19 






2.75±0.75 


232 


5.31 


1500 


12.17 


1.108 


0.712 










26.259±0.19 






4.00±0.50 


54 


5.31 


1500 


12.00 


1.170 


0.662 










25.889±0.24 




Paltani et aL (2007) 


3.50±0.50 


113 


1720 


1700 


3.80 


1.000 


1.000 


-2 91+01* 


-21.49±0.19 


-1.40f 


26.0931016 




26.499l0;09 


Reddv et aL (2008) 


2.20±0.32 


10007 


1925.8 


1700 


3.796 


1.000 


1.000 


-2.76+0-" 

— 0.20 


-20.97±0.23 


-1.84±0.11 


26.439+0-06 

—0.07 


26.256 


27.030 


Sawicki & ThomDson (2006a, b) 


2.20±0.32 


2417 


169 


1700 


5-15 


1.000 


1.000 


-2 52+°-20 


-20.60tO-j8 


^■^"-0.22 


26 320+0 02 


26.314 


26.424±0.03 




2.96±0.26 


1481 


169 


1700 


5-15 


1.000 


1.000 


-2 77+0-" 

^- ' ' -0.09 


-20 go+O'^^ 


-1 43+0-l'^ 
^■^■^-0.09 


26.257±0.01 


26.303 


26.422dz0.03 




4 1 3+0 26 


427 


169 


1700 


5-15 


1.000 


1.000 


-3 07+°-2i 


-21 OO+O-*" 


-1 26+0-'*° 


25 969-1-0 03 


25 965-1-0 03 


26 061 -1-0 07 


ajimasaku et aL (2005) 


5.90±0.30 


12 


767 


1425 


4.358 


1.000 


1.000 










24.44710-11 




Shim et aL (2007) 


3.20±0.14 


1088 


9468 


1550 


3.964 


1.000 


1.000 


-2.81'= 


-20.69° 


-0.83<=f 


26.027 


25.578 


26.067 


Steidel et aL (1999) 


2.96±0.26 


1270 


1046 


1700 


3.796 


1.126 


0.693 


-2.86 


-21.07±0.15 


-1.60±0.13 


26.462 


26.193 


26.711 




4.13±0.26 


207 


828 


1700 


3.796 


1.175 


0.658 


-2.97 


-21.11 


-1.60^ 


26.397 


25.734 


26.640 


Wadadekar et aL (2006) 


2.75±0.75 


125 


5.67 


1850 


4.70 


0.862 


1.249 




-21.14+!!-" 

— 0.15 




26.518 






Yoshida et aL (2006) 


4.00±0.30 


3808 


875 


1500 


3.381 


1.000 


1.000 


-2 84+"-" 
^•°^-0.12 


-1.82±0.09 


26.450 j:E;:°6 


26.343±0.02 


26.968 




4.70±0.30 


539 


875 


1500 


3.381 


1.000 


1.000 


-2 Ql+O " 


-20 


-1.82f 


26.138±0.24 


25.645lH^ 


26.726 


This work 


2.28±0.33 


7093 


857.5 


1700 


3.796 


1.000 


1.000 


-2.25±0.46 
-2.50±0.17 
-2.75±0.21 


-20.50±0.79 
-20.95tO;« 
-21.30±0.35 


-1.05±1.11 
-1.60' 
-1.84^ 


26.519±0.68 
26.601±0.07 
26.633±0.08 


26.276±0.68 
26.293±0.07 
26.305±0.08 


26.601±0.92 
26.864dz0.07 
27.172±0.08 



Note. — Cols. (1) through (5) list the reference, the redshift, the sample size, the area (arcmin^), and the rest-wavelength of measurements (A). Dust extinction correction is provided in Col. 
(6), and corrections to a common cosmology for luminosity and number density are shown in Cols. (7) and (8) . Schechter LF parameters are listed in Cols. (9)-(ll) in units of Mpc~^ for number 
density, and finally the observed luminosity densities (erg s"""" Hz""'" Mpc~'^) are provided in Cols. (12)-(14) for three different limits of integration. 

^ These values were not reported in the paper, but was obtained by fitting their LF with data provided in the paper or by the authors. For lShim et all (]2007l') . the results arc not well constrained 

f 

given the limited number of degrees of freedom. This value was kept fixed, while the other Schechter parameters were fitted. 



